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We study multi-field DBI inflation models with a waterfall phase transition. This transition hap- 
pens for a D3 brane moving in the warped conifold if there is an instability along angular directions. 
The transition converts the angular perturbations into the curvature perturbation. Thanks to this 
, conversion, multi-field models can evade the stringent constraints that strongly disfavour single field 

I ■ ultra-violet DBI inflation models in string theory. We explicitly demonstrate that our model satis- 

fy fies current observational constraints on the spectral index and equilateral non-Gaussianity as well 

£\j _ as the bound on the tensor to scalar ratio imposed in string theory models. In addition we show 

that large local type non-Gaussianity is generated together with equilateral non-Gaussianity in this 
" model. 
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I. INTRODUCTION 



The inflationary scenario has been established as a standard model for the very early universe not only because it 
solves the problems of the standard big bang scenario such as the horizon problem, flatness problem and monopole 
problem, but also because it explains the origin of the almost scale- invariant spectrum of primordial curvature per- 
turbations that seeded the Cosmic Microwave Background (CMB) anisotropies (see, e.g. (lj, l2|). However, there 
are numerous models of inflation that are compatible with the current cosmological observations. Therefore, more 
precise observations, such as those from the PLANCK satellite [3j, will help us distinguish between many possible 
early universe models. 

The origin of innaton, the scalar field that is responsible for inflation, is not specified in many inflation models. 
Dirac-Born- Infold (DBI) inflation 0, 0] motivated by string theory identifies the innaton as scalar fields describing the 
positions of a D-brane in the higher dimensional space in the effective four-dimensional theory. Although DBI inflation 
is well motivated and it predicts interesting features such as the speed limiting effect on the velocity of the scalar fields 
@ and large non-Gaussianity 0, current observations already give strong constrains on the models. In particular, it 
was shown that the Ultra- Violet (UV) DBI inflation models where a D3 brane is moving down the warped throat is 
7-H . already ruled out by current measurements of the tensor to scalar ratio, the spectral index and non-Gaussianity when 
one applies microphysics constraints on the variation of the innaton field in a string theory set-up 

However, in the multi-field DBI models, such stringent constraints can be relaxed. In fact, DBI inflation is naturally 
■ a multi-inflation model as there are six extra dimensions which are the radial direction and five angular directions 
in the internal space. In multi-field models, the trajectory in the field space can have a turn and it converts the 
04 ■ entropy perturbations into the curvature perturbation on superhorizon scales [l3j . It was shown that if there is a 
sufficient transfer from the entropy perturbations to the curvature perturbation, the constraints on DBI inflation can 
. be significantly relaxed However, in order to make definite predictions, it is required to calculate the transfer 

T *~j ' coefficient explicitly with a concrete multi-field potential. 

The potential for the angular directions for a D3 brane in the deformed warped conifold was calculated in Ref. [l8[ 
and the impact of angular motion on DBI inflation has been studied in Ref. fl9| . Ref. [2(| shows that the angular 
directions can become unstable in a particular embedding of D7 brane on the warped conifold and the angular 
instability connects different extreme trajectories. These potentials are calculated assuming that the backreaction of 
the moving brane is negligible and it cannot be applied to DBI inflation directly. However, it is natural to consider 
that a similar transition due to the angular instability happens also in DBI inflation. The potential derived in Ref. [2(| 
has a similar feature to the potential in hybrid inflation. The mass of the entropy field is large initially. As the brane 
moves in the radial direction, the mass becomes lighter. Eventually, it arrives at the point where the entropy field 
becomes tachyonic. Then the inflaton rolls down to the true vacuum along the entropy direction and moves down in 
the radial direction along the true vacuum. 

In this paper, we analyse a two-field DBI model with a potential which has a similar feature as those obtained in 
Ref. [20j |. Using this potential, we will explicitly study predictions for observables such as the spectral index, tensor 
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to scalar ratio and non-Gaussainity and see if we can avoid the stringent constraints that rule out the single field 
UV DBI inflation models. To simplify the calculations, we consider a potential which has an effective single field 
inflationary attractor with a constant sound speed [2l[ before and after the tachyonic instability develops along the 
angular direction. DBI inflation models are known to generate large equilateral type non-Gaussianities Q from the 
bispectrum of the quantum fluctuations of the scalar field before the horizon exit. On the other hand, large local type 
non-Gaussianities [22| (see [23| for a review and references therein) can be generated on super-horizon scales in multi- 
field models because of the conversion from the entropy perturbations to the curvature perturbation. Therefore, 
in general multi-field DBI inflation predicts a combination of the equilateral type and local type non-Gausianities 
and this feature can be used to distinguish DBI models from other inflationary models [24-26]. The presence of 
both equilateral and local type non-Gaussianities in multi-field DBI inflation was first pointed out by Ref. [27j (see 
also Ref. [28]). Ref. [27j considered a model where the conversion of the entropy perturbations into the curvature 
perturbation happens at the end of inflation. In this paper, we consider the case where the conversion happens during 
inflation by a waterfall phase transition due to the instability along the angular direction. 

This paper is organized as follows. In section HH we briefly review the background dynamics of the single field DBI 
model with a constant sound speed. Then, the two-field potential that we consider in this paper is introduced. We 
also show numerical results for the background dynamics with this potential. In section Hill the dynamics of the linear 
cosmological perturbations are studied. We show the results for the power spectrum of the curvature perturbation, 
the tensor to scalar ratio and the equilateral non-Gaussianity in this model and demonstrate that it is possible to 
evade the constraints that rule out the single field UV DBI inflation models. In section ITV1 we compute the local-type 
non-Gaussianity of the curvature perturbation by using the (5N-formalism. In section |Vj we summarise this paper. 
In Appendix A, we review the decomposition in the field space. In appendix B we explain the numerical method to 
analyse the linear perturbations using the decomposition and the <5N-formalism. 



II. BACKGROUND DYNAMICS 



In this section, we first review the background dynamics in the single field DBI inflation model that has a late time 
attractor solution with a constant sound speed. Then, we study the background dynamics in the two-field model with 
a potential that leads to a waterfall phase transition. 



A. The model 



The Lagrangian for the multi-field DBI inflation is given by 

P{X IJ , cp 1 ) = P(X, <ff) = -j^ yi-2fW)X -lj-V (cp 1 ) , (1) 

where ch 1 are the scalar fields (I = 1,2,...), f(4> 1 ) and V(4> 1 ) are functions of the scalar fields and X is defined in 
terms of the determinant 

V = det(£5 - 2/X/) 

= i - 2fGuX IJ + 4f 2 x l /xj ] - sfxfxjx^ + mfxfxjxZx? , (2) 



as 



where 



X IJ = -\d^d^ J , (4) 

and G IJ is the metric in the field space. Note that /(0 J ) is defined by the warp factor /i(^ / ) and the brane tension 
T 3 as 
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The sound speed is denned as 



P 



\ Px + ^P,xx 



= \^l-2fX, 



(6) 



where means the partial derivative with respect to X. Note that X coincides with X = GijX in the homogeneous 
background because all the spatial derivatives vanish. From the action (JXJ) , we can show that 



P) 



In this paper, we consider the Einstein-Hilbert action for gravity and hence all equations of motion are derived from 
the action 



S=^Jd 4 x ^R + 2P(X IJ ,cj ) I ) 
where we set 8ttG — 1 and is the four dimensional Ricci curvature. 



(8) 



B. DBI inflation with a constant sound speed 

Let us consider single field DBI inflation with the potential Vi(<f>) and warp factor /(</>) where cj> is the scalar field. 
Then, the field equation is given by 

jS + 3 Hj> - ^ + esV^ - (1 ~ Cs)2 4# =0- (9) 

In [2l|, it was shown that when Vi(<fi) and /(</>) are given by 

V!(ct>) = Vo<p- q , (10) 

/(^) = M q+2 , (ii) 

with constants Vo, fo and q, Eq. (|9]) has a late time attractor inflationary solution with a constant sound speed that 
is given by 



IGfoVo + 3 



(12) 



Throughout this paper, as a concrete example, we consider the case with q = 4. This attractor solution is potential 
dominated, which means that the potential term is much larger than the kinetic term along the attractor solution. If 
we define the slow-roll parameters as 

€ = -JT^ V= lTe> S= Hc~ s - (13) 

s vanishes for this attractor solution and also e and rj are much less than unity because it is a potential dominated 
attractor solution. 



C. Two-field model 



We investigate the following two field model. We define the field as a radial direction and define the field \ as an 
angular direction in the warped throat. The field space metric is then given by 



Gu = AiSu, 



(14) 
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where = 1, A x = <f> 2 and 8u is the kronecker delta. With this field space metric, X becomes 

x = \U 2 + 4?x 2 ), (15) 



2 

in the homogeneous background. The Friedmann equation is given by 

1 ( 1 



31,1 'TwAi- 1 )^- (16) 

By varying the action (j8]) with respect to the fields, we obtain the equations of motion for the fields as 

+ 3 jty _ ^ _ fa 2 + CsVa> _ Q^LtA = 0, (17) 



X + 3Hx--x)+2<Ma + c.V, x = 0. (18) 



If we specify the potential and the warp factor, the background dynamics are determined by solving Eqs. (|17[) and 

mi. 

In DBI inflation, the scalar fields describe the positions of a brane in the bulk. The explicit form of the multi-field 
potential depends on the details of the geometry of the warped conifold and various effects from the stabilisation of 
moduli fields. In Ref. [20(, an example of the multi-field potential that has a similar feature with hybrid inflation 
was obtained. In this potential, the inflaton rolls down along the radial direction first. Eventually, it arrives at the 
transition point where the entropy field becomes tachyonic. Then the inflaton rolls down to the true vacuum along 
the entropy direction and moves down in the radial direction along the true vacuum. 

In this paper, we investigate a two field potential which captures the essential feature of the potential derived in 
string theory as described above. We assume the radial field has the form of the potential with a constant sound 
speed as is discussed in the subsection III Bl to simplify the calculation. The two field potential is given by 

v{^ x ) = \Kx 2 -xlf + g[fj ' + ^, (19) 

Let us assume that the inflaton starts rolling down in the radial direction (j) with a small deviation from x = 0. Then, 
in the early stage when x <C 1, the potential is effectively a single field potential for </>; 



V{<t>,x)~\\xi + % (20) 



At this stage, if we assume 



the effective potential becomes 



2 



A^«^, (21) 



V cS ,i~^- (22) 

Because this is in the same form as the potential (fTU|). there is a late-time attractor solution with a constant sound 
speed given by Eq. (|12j) . As the inflaton rolls down in the radial direction, there appears a waterfall phase transition 
as is shown in Fig. Q] where the inflaton rolls down to the true minimum of the potential (|19[) . This can be seen more 
clearly by rewriting the potential (fl9l) as 



v(4>,x) = \\ 



.2 I „2 _9_ 



X - \ Xo ,j. 2 



"I 2 / „2 



In this form, x appears only in the first term. We can clearly see that x = is the minimum in the x direction when 
4> 2 < <?/(Axo)j while x 2 — Xo ~~ 9/{^4> 2 ) becomes the minimum in the x direction when (f> 2 > g/(XXo)- Therefore, 



5 




Log(V) 




FIG. 1. The potential given by Eq. (|19[). In order to show the feature of the potential clearly, we plot log V(<j>, x)- The 
parameters are chosen as A = 1.2 x 10~ 6 , %o = 0.004, Vo — 5 x 10 -12 and g = 3 x 10~ 9 . For small (f> there is a minimum at 
X — but for large (j> t rne vacua appear and the field rolls down the to the true vacua causing the waterfall phase transition. 



<f> 2 = gj (Axo) is the critical transition value for <p. The effective potential in the true vacuum with \ 2 — Xo ~ 9/(^4> 2 ) 
is given by 

2 



Thus if we assume 



V(<j>,x)~V 1- 



g{xo/4>Y 



9 \ 1 ( Xo 



2XV 



-Z+5 — 



gxl 4> 2 



V (1 - g 2 /2XV Q ) 1/cb 4 Vo l-g 2 /2XV Q 
the effective potential in the true vacuum becomes 

9 2 \ 1 



«1, 



Kff,2 ~ V l - 



2AV 



Again, this is in the form ~ <fi 4 and we have a late-time attractor with a constant sound speed 



+ 3 



where 



V = V 1 



2XVn 



(24) 



(25) 



(26) 



(27) 



(28) 



Now we show our numerical results for the background dynamics. We choose the parameters as follows; A = 
1.2 x 1(T 6 , xo = 0.004, V = 5 x 1CT 12 and g = 3 x 1CT 9 . The warp factor is given by Eq.([TI]) with / = 1.2 x 10 15 . 
There parameters are chosen so that all the observables satisfy the current observational constraints. 

Firstly, the left panel of Fig. [2] shows the dynamics of the inflaton in the x direction. Before the transition happens, 
the potential has its minimum at x — in the x direction. Therefore, regardless of the initial conditions, the inflaton 
rolls down the potential and the value of x approaches unless the transition occurs while x is still large. We use the 
c-folding number N — In a as time. We normalise the e-folding number so that the transition finishes after N = 60. 
In this example, x 1S sufficiently small at N = 10 and we have an effective single field dynamics until around N ~ 15. 
As we expected, the inflaton rolls down to the true vacuum in the transition, which occurs during 20 < N < 55. After 
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the transition ends, it rolls down along the true vacuum. Notice that the true vacuum is not along a constant \ nne 
but the value of x along the true vacuum Xtrue is a function of (f>; 



Xtr 



xl 
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(29) 



which can be obtained from Eq. (|23|) . However, as we see in the right panel of Fig. [2j the trajectory along the 
true vacuum curves slowly so that the coupling between the adiabatic mode and the entropic mode can be ignored. 
Actually, using the value of the Hubble parameter in this model which is H ~ 2.8 x 10~ 8 we can roughly estimate how 
many e-folds we need after sound horizon exit of the mode which we are considering [29j assuming instant reheating. 
In our model, it is around 60 e-folds. We can see that the transition ends within 60 e-folds after sound horizon exit if 
we consider modes that exit the sound horizon in the effective single field regime (10 < N < 15). We assume inflation 
ends after the transition by some mechanisms such as an annihilation of D brane with anti-D brane. 

The right panel of Fig. [2] shows the sound speed. Before the transition, the sound speed slowly changes. This is 
because the condition (|2"Tj) is not fully satisfied. On the other hand, after the transition, we can clearly see that the 
sound speed is almost constant. The sound speed changes the most during the transition and the slow-roll parameter 
s takes the largest value during the transition. However, the largest value of s is still around —7 x 10~ 4 . Actually, 
as is shown in Fig. [31 all the slow-roll parameters are always much smaller than unity even during the transition. 
Therefore, the slow-roll approximation always holds in this model and the sound speed is almost constant even during 
the transition. 
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FIG. 2. Left: The dynamics of the x field. In the early stage (10 < N ^ 15), the inflaton rolls down the potential almost in 
the 4> direction with x ~ 0. The transition occurs during 20 < N < 55. After N ~ 55, the inflaton rolls along the true vacuum. 
Right: The sound speed. The sound speed is almost constant because it changes very slowly even during the transition. 



III. LINEAR PERTURBATION 



In this section, we study linear perturbations in the two field DBI model introduced in section |n] We first derive 
the coupled equations for adiabatic and entropy perturbations. Secondly, we present the numerical results for the 
power spectrum of the curvature perturbation. Then, we show that the constraint on the tensor to scalar ratio in 
the single field DBI inflation is incompatible with current observations and show how this constraint can be relaxed 
in the multi-field models. Finally, we demonstrate that the model considered in this paper can actually evade the 
constraints. 
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FIG. 3. Left: The slow-roll parameter e. It is always smaller than 1 x 10 -3 even during the transition. Middle: The slow-roll 
parameter r/. Its absolute value is always smaller than 2 x 10 -3 . Right: The slow-roll parameter s. \s\ takes the largest value 
—7 x 10~ 4 during the transition. However, its absolute value is still much less than one. 



A. Field perturbations 



We introduce the linear perturbations of the fields and Q x as 

<t> = MQ + Q+(?,t) , x = xo(t) + Q x {x,t) 



(30) 



As discussed in appendix [A] we decompose these perturbations into the instantaneous adiabatic and entropy pertur- 
bations as 



Qf = Qaet + Q s ef, 

where e^. and are the adiabatic and entropic basis, respectively. They are defined as 

h 6 1 <l>x 



S 2 + cf> 2 x 2 



> 2 + (j) 2 X 2 



(31) 
(32) 

(33) 



e x 



1,2 _|_ 02^2 



i 2 + 4> 2 x 2 



(34) 



For the analysis of perturbations, it is convenient to use the conformal time r = J dt/a(t) and define the canonically 
normalized fields as 



V a = —Qcr, V s = —Q s 

C s C s 



The equations of motion for v a and v s are obtained as 



(*£)' 



Vs = 0, 



(35) 



(36) 



< + iv' a + [ c 2 s k 2 - — + a 2 fx 2 ) v s - -fr a = 0, 



where the prime denotes the derivative with respect to r and 



i + c 2 s )p a -4& 2 P <Jts 



(37) 



(38) 
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1 ~ 9 

^S 1 ,SS . n 1 t 



(39) 



with 



S2X, P s = PjeiV^s, P^xs = PxAV^s, P ss = (ViVjP) e{ 



e J s c a , 



(40) 



(41) 



where T>j denotes the covariant derivative with respect to the field space metric Gjj. We numerically solve Eqs. (j3"6")l 
and (|37p to compute the curvature perturbation. 



B. Curvature perturbation 

In this subsection, we show how we set the initial conditions for Eqs. (|3f)|) and (j3"7| to calculate the power spectrum 
of the curvature perturbation. If the trajectory is not curved significantly, the coupling £/ 'aH becomes much smaller 
than one. In Fig.[H we see that £,/aH is still much smaller than one before the transition starts around N ~ 20. Also, 
the slow-roll conditions are satisfied as is shown in section III CI and this means that H, a and c s change very slowly 
with time compared to the Hubble scale so that the approximations z" j z ~2/t 2 and a" /a ~ 2/r 2 hold. Therefore, 
we can approximate the Eqs. (|3"6"f and (|3"T1) as Bessel differential equations before N ~ 20 (see appendix [B] for the 
entropy perturbation). Then, the solutions with the Bunch-Davis vacuum initial conditions are given by 



y/2kc s 



-e ~ lkCsT [ 1 



(42) 



1 



Vsk 



__ e -ikc s r , 1 



\J2kc s V kc s r / 

when (J,g/H is negligible for the entropy mode. Then, the power spectra of Q„ and Q s are obtained as 

H 2 „ H 2 



Qa 



Att 2 c s ' * Qs ~ 4ir 2 c s ' 

which are evaluated at sound horizon crossing. The power spectrum of 1Z at sound horizon crossing is given by 



27T 2 ' 



Pn, = -^\n\ 



k 3 M 2 


H 4 




H 2 


2n 2 z 2 


Att 2 ct 2 

* 


* 


8ir 2 ec s 



(43) 



(44) 



(45) 



where the subscript * indicates that the corresponding quantity is evaluated at sound horizon crossing kc s = aH. 

We investigate a mode which exits the sound horizon at N ~ 10 where we have an effectively single field dynamics. 
As stated above, the coupling £/aH is much smaller than unity around sound horizon exit. Also, as we can see in 
Fig. 21 \{j, 2 /H 2 \ is also much smaller than unity around sound horizon exit at N ~ 10. The mass also changes very 
slowly. If we define a quantity 



Mr. 



(46) 



which quantifies how rapidly the mass of the entropy perturbation changes, we can see in Fig. 2] that M c is smaller 
than unity for at least 5 e-folds after sound horizon exit. Therefore, we can set the initial conditions for Eqs. (|36|) 
and (l37l) by the solutions (|42|) and (|43l) . Note that we set the initial conditions at N ~ 7 when the mode which we 
consider is still well within the sound horizon. 

We treat v a and v s as two independent stochastic variables for the modes well inside the sound horizon as in [30l |. 
This means that we perform two numerical computations to obtain Pn for example. One computation corresponds 
to the Bunch Davis vacuum state for v a and v s is set to be zero to obtain the solution 1Z\. Another computation 
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FIG. 4. Left: The coupling ^/aH. Before and after the transition the coupling is small but it becomes large during the 
transition converting the entropy perturbation into the curvature perturbation. Middle: The mass of the entropy perturbation 
divided by the Hubble parameter squared. Its absolute value is much smaller than one around sound horizon exit at N ~ 10. 
Right: The quantity M c that shows how the mass of the entropy perturbation changes. It is also smaller than one around 
sound horizon exit. 



corresponds to the Bunch Davies vacuum state for v s and v a is set to be zero, in which case we obtain the solution 
1Z2- Then, the curvature power spectrum can be expressed as a sum of two solutions; 

Pn = ^(\n 1 \ 2 + \n 2 f). (47) 

This procedure is applied to all the numerical computations in this paper. 

Fig. [5] shows the numerical results for the power spectrum of the curvature perturbation. As we can see from the 
small value of e, the Hubble parameter H changes very slowly. It changes only by a few percent between TV ~ 10 and 
N ~ 70. Substituting H ~ 2.8 X 10~ 8 , c s ~ 6.2 x 10~ 3 and e = 6.4 x 10~ 4 into Eq. @5J), we obtain the value of the 
curvature power spectrum after sound horizon exit as 



P - H 



2.5 x 10~ 12 . (48) 



This should coincide with the result of the numerical solution. Actually, in the left figure of Fig. [5j it is shown that 
Pji becomes almost constant soon after sound horizon exit and it is given by Pji ~ 2.5 x 10~ 12 . It does not change 
significantly until N ~ 15. When the trajectory in field space curves, the curvature perturbation is sourced by the 
entropy perturbation and it is enhanced. We can actually see that the power spectrum of the curvature perturbation 
is enhanced by a factor of ~ 9 x 10 2 during the transition in the right figure of Fig. [5] After the transition, it takes a 
constant value Pji ~ 2.3 x 10 -9 , which is compatible with the CMB observation [l(. 
If we express the final curvature power spectrum as [16| 



n 



cos 2 9 ' 



(49) 



the enhancement is quantified by the function cos 2 8. In this model, we have cos 2 6 ~ 1.1 x 10 3 . When \/j, 2 /H 2 \ is 
much smaller than unity, from Eq. (|49p . the spectral index n s is given by 

n s — 1 = ^1"^ = —26* - 77* — a* sin 29 — 2/3, sin 2 0, (50) 
a In A; 

where 

-SP'-H-jMrf + D^-T- <51) 

where only the leading order terms in the slow-rolling approximation are kept in the expression for f3. Also, the 
parameter which quantifies the equilateral non-Gaussianity is expressed as 



iNL - 108c 2CO& ° 3c? • (52) 
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FIG. 5. Left: The power spectrum of the curvature perturbation becomes constant a few e-folds after sound horizon exit which 
is around -/V ~ 17. Right: We can see that the power spectrum of the curvature perturbation is enhanced during the transition. 
Then, it becomes constant again after the transition. Note that both of these figures describe the contribution from the real 
part of v a h and the contribution from the imaginary part shows the same behaviour. Therefore, it is sufficient to show only 
these figures in order to know its behaviour because the power spectrum of the curvature perturbation is just a sum of those 
two contributions. 



C. Gravitational waves constraints 

In Ref. it is shown that the single field ultra-violet (UV) DBI inflation is disfavoured by observation as follows. 
Firstly, the tensor to scalar ratio is related to the field variation Acf> by the Lyth bound 



M 2 P \ANJ 



(55) 



where N = J dtH. Because the field variation A</> corresponds to the radial size of the extra dimensions, it is 
constrained by the size of the extra dimensions. From Eq. (|55p . the upper bound on A(f> gives the model independent 
upper bound on the tensor to scalar ratio for standard UV DBI inflation. The bound is typically given by 

r < 10" 7 , (56) 

when we assume the minimum number of e-foldings that could be probed by observation as AN ~ 1. Secondly, the 
lower bound on the tensor to scalar ratio in the single field UV DBI inflation is derived in the following way. The 
relation (l50l) can be rewritten as 



3i/srv / 



4cos 3 9 Hf 



a*sin29 + 2/3*sirr 6, (57) 
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by using the Eqs. (IT3l) and (|52j) as shown in Ref. [l6|. Note that a term proportional to c 2 s» is neglected because 
both c s and s* are small. For the single field UV DBI inflation, we have / > and = 0. This gives 

4 

r > — ^=^= (1 - n s ) (single field), (58) 

from Eq. (|58|) . The amplitude of the equilateral non-Gaussianity is constrained as 

- 151 < fiff 1 < 253 at 95% C.L., (59) 

from WMAP5 and the best-fit value for the specrtral index is n s ~ 0.982 From those values, we can obtain the 
lower bound on the tensor to scalar ratio as 

r > 10~ 3 . (60) 

Clearly, the lower bound (|55|) is not compatible with the upper bound (|60p . This is why single field UV DBI inflation 
is disfavoured by observation. 

These constraints are relaxed when we consider the multi-field models. The upper bound is relaxed because we 
have angular directions and the field variation is not only determined by the radial coordinate. More importantly, 
the lower bound is relaxed significantly. In Eq. (|57[) . the last two terms become important if there is a transfer from 
entropy to adiabatic modes (6 ^ 0). In the model considered here, the curvature perturbation originated from the 
entropy perturbation dominates the final curvature perturbation and the spectral index is indeed determined by the 
mass of the entropy mode 



2^ 



3H 2 



1 - 0.972, (61) 



which is compatible with the WMAP observation. This value has also been confirmed in our numerical computations 
for the curvature perturbation. Thus there is no longer the lower bound for the tensor to scalar ratio. The tensor to 
scalar ratio is obtained by substituting e ~ 6.4 x 10~ 4 , c s ~ 6.2 x 10~ 3 and cos 2 9~l.lx 10~ 3 into Eq. ([5i| as 



16ec s |„cos^e ~ 7.0 x 10~ 8 . (62) 



This is compatible with the upper bound (|56l) . 



IV. NON-GAUSSIANITIES 

In this section, we calculate non-Gaussianities of the curvature perturbation. In addition to the equilateral type 
non-Gaussianity, the transition may generate local type non-Gaussianity. The local type non-Gaussianity can be easily 
calculated by the SN formalism. We first briefly review the <5N-formalism [3l|, |32| and compute the power spectrum 
of the curvature perturbation in order to confirm the accuracy of the SN formalism in our model. Then, we use it to 
compute the non-Gaussianities of the primordial curvature perturbation. 

A. 5N formalism 

In the (5N-formalism, the curvature perturbation on the uniform density hypersurface £ evaluated at some time 
t — tj is identified as the difference between the number of e-folds iV^ (i/ , t j , x) and N (tf,ti) where iV^ (£/ , t j , x) is 
the number of e-folds from an initial flat slice at t = t{ to a final uniform density slice at t = tt and N Q (tf,ti) is the 
number of e-folds from an initial flat slice at t = ti to a final flat slice at t = tf\ 

((t f ,x.)~dN = N A (t f ,t h x)-N (t f ,ti). (63) 

We can ignore the dependence of N on the time derivatives of the fields when the slow-roll approximations hold 
because the equations of motion for perturbations are reduced to first order differential equations. Then, we can 
expand SN in terms of the field values at the sound horizon crossing up to the second order 



SN [ <Pa,^^s, <Ps)^SN ((f> a ,fa) (64) 

- N,^Q a + N^Q S + \N KK Ql + \n M ,QI > (65) 
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where <\> a denotes the scalar field in the instantaneous adiabatic direction which is almost equivalent to <f> before the 
transition and <f> s denotes the scalar field in the instantaneous entropic direction which is almost equivalent to x before 
the transition (see Fig. [2]). Note that we ignore the cross term Q a Q s because the two fields are independent quantum 
fields around sound horizon exit before the transition. 

We compute the power spectrum of the curvature perturbation at t = tf after the transition by the SN- formalism 
taking U to be well before the transition when the curvature power spectrum is constant after sound horizon exit. In 
this case, as is shown in Fig.[l)J the transition occurs at different <f>„ if we perturb the initial field values in the entropic 
direction. Because taking perturbations in the initial field values changes the trajectory unlike single field cases, we 
need to be careful about the definition of the final slice. Let us define as the field values along the unperturbed 
trajectory (solid line) and as the field values along the perturbed trajectory (dotted line) in Fig. [6] We define SN 
as 

SN = J Hdi~ J Hdt, (66) 

where the tilde denotes the quantities with the perturbed initial conditions. For example, if we perturb the initial 
field values to the entropic direction, the perturbed initial field values (j> a \ti) corresponds to the position in the field 
space as 

4 > i(t i )=<f>i(t i )+Q s el, (67) 

with el given in Eqs. (|3"3")) and (|34p. Note that if in Eq. (f6"6"|) is the time when takes the same value as 4>i(tf) 
well after the transition. Because both trajectories merge into the attractor solution in the true vacuum after the 
transition as is shown in Fig. [6j we take the final slice so that both unperturbed and perturbed trajectories have 
the same field values <j)f and \f on the final slices. Because we can determine all the phase space variables if we 
know the values of the fields in the slow-roll case, this means that we have the same cj), 4>, \ and \ on the final slices 
which results in the same H (i.e. uniform density) from Eq. (|16[) . By using the definition of 5N in Eq. (|66p. we can 
numerically compute the quantity 

N, s = ^ , (68) 

where we make Q a sufficiently small so that the value of N t< j, a does not depend on the value of Q a . Because the 
contribution from N t< j, a in Eq. (165[) is negligible in our model, we obtain 

P n = P c ~N%P Qs \ t=u , (69) 

where Pq 3 is given by Eq. (|44[) . Note that the curvature perturbation on the uniform density hypersurface coincides 
with the comoving curvature perturbation on super-horizon scales (33j . The numerical result shows that the <5N- 
formalism successfully predicts the value of the final curvature perturbation P-n ~ 2.3 x 10 -9 within a few percent 
error in this model. 



B. Non-Gaussianities 



Now, we can compute the non-Gaussianities of the curvature perturbation using the SN formalism as follows. The 
bispectrum of the curvature perturbation is defined as 

(C(£)C(£)C(£)) = (2TT) 2 B c (k 1 ,k 2 ,k 3 )S (3 \k 1 +k 2 + k 3 ). (70) 

We can rewrite Eq. (1701) as 

(C(*i)C(fe)C(fe)) - N^N% (Q CT (fci)Q s (fcl)Q s (fcl) + (perm)) 

+ ^N%N Mm (Q s {ki)Q s {k 2 ) (Q s *Qs) (fca)) + two perms. (71) 

Here we used the fact that N^^ is much smaller than A^ s in Eq. and the dominant bispectrum of the fields 
is coming from the mixed adiabatic and entropy contributions. Note that the star * denotes the convolution and 
correlators higher than the four-point were neglected in the above equation. From Eqs. (fTDj) and (|7ip . we obtain 



B c {k l7 k 2 ,k 3 ) = —L^N^BQ^k^h^h) 



' 4 ^TT^a^7' (72) 



13 




FIG. 6. Left: Two trajectories in the field space obtained by perturbing the initial value of the entropy field. The solid line 
describes the unperturbed trajectory while the dotted line describes the perturbed trajectory. Right: The transition occurs at 
different values of (f). Again, the solid line is the unperturbed trajectory while the dotted line is the perturbed trajectory. 



by using Eq. (|(r9"|) where the bispectrum of the scalar field perturbation BQ a (ki, k2,k 3 ) is denned as 

[Q*(ki)QAk2)Qv(k 3 )) = (2^) 2 B QCT (fc 1 ,fc 2 ,fc 3 )^(fc; + k 2 + k 3 ). (73) 

Here we used the fact that the symmetrised mixed bispectrum has the same shape as the pure adiabatic bispectrum 
[13 ] and N t rt >4> /N. l p s ~ cos 9 <C 1. The non-linear parameter /nl is define as (34) 



_ 5 1 1 /• v; B C 

= 



If the non-Gaussianity is local, we can express the curvature perturbation as 

C - Cn + |/ N L (C ~ (O) , 

where obeys Gaussian statistics. From Eqs. ([72)1 and (j?4"|) . the nonlinear parameter /nl is expressed as 

5 Na.a. 



£ 1 r\ requil 

/nl = cos z QfJr 



i requil 

where f* L 



e = o 6 N 2 ^ ' 

is the equilateral non-Gaussianity parameter in the single field model defined as 



(74) 



(75) 



(76) 



equil 
NL 



5 nA 3 b, 



Qa 



e=o 6 Sjfc 3 3 ^N^P^ 



35 1 
108 cf' 



The first term on the right hand side of Eq. (|76l) , 



requil _ 2 q r equil 
J NL — Lob ° JNL 



e=o 



(77) 



(78) 



comes from the bispectrum of the quantum fluctuation of the scalar fields generated under horizon scales and gives 
rise to the equilateral non-Gaussianity, Eq. ([5"2]). We can see that small cos0 suppresses f^if- ^ n our m °del, the 
equilateral non-Gaussianity is obtained as 



requil 
JNL 



-9.5. 



The second term on the right hand side of Eq. (I7B1 , 



f local 
INL 



5 

6 N 2 



(79) 



(80) 
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is generated even if the field perturbations at the horizon crossing are Gaussian. This is called local type non- 
Gaussianity. Numerically, we can compute the second derivative of N with respect to cf> 8 as 

N (#(fr) + Q s ej) - 2N + N - Q s ej) 

The result of our numerical computation shows that we have 

&T' = !^^~40.1 (82) 

in our model (see appendix IBl for the details of the numerical computations). Note that /j^£ a/ becomes constant after 
the transition. This is compatible with the current WMAP observation [l[ 

- 10 < < 74 at 95% C.L. (83) 



V. SUMMARY AND DISCUSSIONS 



DBI inflation is the most plausible model that generates large equilateral non-Gaussianity. However, single field UV 
DBI models in string theory are strongly disfavoured by the current observations of the spectral index and equilateral 
non-Gaussianity It has been shown that these constraints are significantly relaxed in multi-field models if there is a 
large conversion of the entropy perturbations into the curvature perturbation. In this paper, for the first time, we 
quantified this conversion during inflation in a model with a potential where a waterfall phase transition connects two 
different radial trajectories (see Fig. [1}. This type of potential appears in string theory where the angular directions 
become unstable in the warped conifold, which connects two extreme trajectories [2(|. We demonstrated that all the 
observational constraints can be satisfied while obeying the bound on the tensor to scalar ratio imposed in string 
theory models. The large conversion also creates large local type non-Gaussianity in general. In our model, this is 
indeed the case and we expect that large equilateral non-Gaussianity is generally accompanied by large local non- 
Gaussainity in multi-field DBI model. This prediction can be tested precisely by upcoming data from the Planck 
satellite. 

There are a number of extensions of our study. In this paper, we study a toy two-field model. It would be important 
to study directly the potentials obtained in string theory to confirm our results although it is a challenge to compute 
the potential when the sound speed is small. The curve in the trajectory in field space, which is essential to evade 
the strong constraints in string theory and responsible for large local non-Gaussianity, can be caused not only by the 
potential but also by non-trivial sound speeds [28j . This happens in a model with more than one throat for example 
where there appear multiple different sound speeds. It would be interesting to compare the two cases to see if one 
can distinguish between them observationally. Finally, it has been shown that the multi-field effects enhance the 
equilateral type trispectrum for a given j^?™' [36[ and its shape has been studied in detail 37]. Moreover, it was 

shown that there appears a particular momentum dependent component whose amplitude is given by /^£ a '/^j 
Thus the trispectrum will provide further tests of multi-field DBI inflation models. 
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Appendix A: Adiabatic and entropy perturbations 

We summarise how the adiabatic and entropic bases arc defined here. We consider the linear perturbations of the 
scalar fields defined as 

cj>\x,t)=^{t) + Q I {x,t). (Al) 
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We decompose the perturbations into the instantaneous adiabatic and entropy perturbations where the adiabatic 
direction corresponds to the direction of the background field's evolution while the entropy directions are orthogonal 
to this For this purpose, we introduce an orthogonal basis e„(n = 1,2. ..AT) in field space. The orthonormal 

condition in general multi-field inflation is given by 

P.xueiei = 8 nm , (A2) 

so that the gradient term P 'x IJ diQ d l Q J is diagonalised when we use this basis. Here we assume that P x IJ is 
invertible and it can be used as a metric in the field space. The adiabatic vector is 



e[ = : , (A3) 

^P X JK<t>-><t> K 



which satisfies the normalization given by Eq. (|A2j) . The field perturbations are decomposed in this basis as 

Q 1 = Quel (A4) 

For multi-field DBI inflation, using the relation 

• • • • 1 — r 2 .... 2X 

P X ;,0V = c s G/j0V + — ^G/kGjl<^W = — (A5) 

IX c s c s 

we can show that the adiabatic vector is given by 



(A6) 



This implies that 

Gueiel = c s , (A7) 



P x u = e,Gjj + l—^-G IK G.j L e? ef . (A8) 



Substituting Eq. (|A8Jl into Eq. (|A2[) . with m = 1, we obtain 

G IJ e[e J n = c s 5 nl . (A9) 
Substituting Eqs. (|A8]l and (|A9|) into Eq. (|A2|) . we obtain 



1 1 - c 2 

G/je^e^ = —5 mn -<W#ni- (A10) 

Cs C s 

For two-field models with Gu (4> K ) — Ai (4> K ) Su where A<p — 1 and A x — </> 2 , from Eq. (|A6|) . the adiabatic vector 
is obtained as 



(All) 

/ 2 + </) 2 X 2 d<j> 2 + 4> 2 x 2 , 
From the orthogonal condition (|A10j) . the entropy vector satisfies 

Gueiei = etc* + 2 e*e* = 0, (A12) 



G /je ^ = (ef) 2 +0 2 (ej) 



(A13) 



which leads to 
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Appendix B: Numerical method 



In this section, we explain how the JN-formalism is used in the numerical computations. In the single field case, we 
just need to perturb the initial conditions along the trajectory. In the numerical computations, it is easy to compute 
4> a (ti + St) because the trajectory with the perturbed initial conditions is the same as the one with the original initial 
conditions. 

However, in the two-field case, the trajectory with the initial conditions peruturbed in the entropic direction is 
different from the one with the original initial conditions as is shown in Fig. |6] Although the perturbed initial values 
of the fields are defined in Eq. (|57]). once we set the value of Q s , we also need to know how to perturb the values of 
the time derivatives of the fields. This is because we need to set the values of all the phase space variables in order 
to solve the second order differential equations numerically. From Eq. (|67j) . we obtain 



fc{t i ) = K{U) + Q s e I s +Q s e I s . (Bl) 
From Eqs. and the derivatives of the entropy basis vector are obtained as 

(B2) 
and 




f <<t> (j> sH a <p p(t) 



(B3) 



3 3 \ 2/ 2 & <f> p(t) J ' 
Here new variables p(t) and q{t) are defined as 

p(t) = -V7^, q(t) = -y/f4>X, (B4) 

so that the sound speed is expressed as 



c s = Vl-p(t) 2 -q(t) 2 - (B5) 



Given that we obtain the numerical values of ej and using Eqs. (|F32|) and (|B3I) . we now know how to perturb all 
the phase space variables (</>, <j), x, x) from Eq. (|Blj) if we know the value of Q s . Because we set the value of Q s , we 
can determine the value of Q s if there is a relation between Q s and Q s . Actually, before the transition, it is possible 
to obtain the analytic solution for v s and hence the solution for Q s from Eq. As mentioned in section Hill the 

approximations z" / z ~ 2/t 2 and a" /a ~ 2/r 2 hold because the slow-roll approximation holds and the coupling £ is 
negligible before the transition. Therefore, Eq. ((57)) is approximated as 



which can be rewritten as 



where 



v" + (c 2 s k 2 - 2/t 2 + a 2 fi 2 s ) v s ~ 0, (B6) 
d 2 v s ldv s ( 9/4- p? s /H 2 



df 2 f df 



0, (B7) 



v s = JL , f = -c s kr. (B8) 

Note that we regard c s as a constant because the slow-roll parameter s is much smaller than unity as we showed in 

Sec. Mil Then, if we approximate fi s to be a constant, we have the analytic solution for Eq. (|B7[> because it is the 
Bessel differential equation. By choosing the Bunch-Davies vacuum initial condition, we obtain 

v s = ^ e *C".+V2W2 (_ r ) V2 H W hCskT) . ( B9) 

where 

us = ^9/4-fi 2 /H 2 , (BIO) 
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and H^} is the Hankel function of the first kind. In the super-horizon limit —c s kr << 1, using the asymptotic form 
of the Hankel function 

H£ , ( - , tT) _ i £M(_i_)"-, (B11) 

we have 



-W7T 



; *K+i/2W2 / _q/2 T (vs) f 2 V 
7r \~c s kT J 



(B12) 



and 



Q s = -^ ^Wi ( _ r) 3/2 £M ( 2 \ ^ 
2 7T \ —C S KT J 

where we used the relation r ~ —1/aH during slow- roll inflation. From Eq. (|B13|) . we obtain the first derivative of 
Q s in terms of Q s as 



3 

s(l - v s ) - e + v s - - 



HQ S . (B14) 



Let us show some numerical results in the model introduced in section III CI As is shown in Fig. [21 the transition 
begins around N ~ 20. We take the initial hypersurface around N ~ 12 where we can evaluate Pq^ with Eq. (|4"4")l 
because the trajectory is still effectively a single field one after sound horizon crossing. On the initial hypersurface, 
we set the values of all the phase space variables as (j){ti) ~ 7.09, cj>(ti) ~ 8.11 x 10~ n , X (ti) — 5.06 x 10~ 5 ,x(£i) — 
3.42 x 10~ 14 . If we set S X = 10~ 7 , we obtain 

^ = ^TIT^ = -2-12 x 10- 9 , (B15) 

from Eq. (|67p where ef(U) and e*(ti) are obtained numerically. Then, we also have 

in- 7 

Q s = = -5.59 x 10" 8 . (B16) 

ef{U) 

Using Eq. (|B14|) . we obtain Q s ~ -4.29 x 10~ 17 . Actually, we can also obtain the first derivative numerically 

— + H(s - 1) ) Q s ^ -3.07 x 10~ 17 , (B17) 

V s J 

from Eq. (1351) using the numerical values of v s and v s ■ We can see that the values of Q s obtained in both ways are 
the same with around 40 percent error. This error comes from the fact that the solution Eq. (|B14p is exact only if fi s 
is perfectly constant and the slow-roll parameters are zero. Using Eqs. (|B1[) . (|B16|) . (|B17[) and the numerical values 
of el and ef obtained by using Eqs. ([33)1 . (jM)) . (|B2|) and (|B3|) . we obtain the first derivatives of the fields as 

8j> = i a (U) - 4> a (ti) ~ -1.47 x 10~ 17 , (B18) 

and 

6 X = xAU) - Xa{h) * 5.35 x 10- 17 . (B19) 

Using Sx — 10 -7 , Eqs. (|B15I) . (|B18I) and (IB19[) . we can now perturb the initial conditions. Using these initial 
condition, the first derivative of TV with respect to <j) s is obtained using Eqs. (jBfJJ) and (|6"8"]l as 

N t<t>s ~ 7.82 x 10 2 . (B20) 

Then, from Eq. (I69[) . the power spectrum of the curvature perturbation is given by 

Pk = Pc- N\P Q , \t= u * 2.29 x 10" 9 , (B21) 
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where we used the numerical result for the power spectrum of the entropy perturbation 

P r 2 \v , I 2 

P Qa \ t=u = f^^L ~ 3.74 x 10- 15 . (B22) 

This coincides with the value Pr ~ 2.3 x 10 -9 with less than one percent error, which is obtained directly by solving 
the equations of motion for the linear perturbations numerically. This confirms the validity of the <5N formalism in 
this model. In a similar way, we perturb N and obtain the second order derivative using Eq. (|8Tj) . Then it is possible 
to compute the local type non-Gaussianity. 
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